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Abstract 

We prove that a flow of cold collisionless particles from all directions in and out of a region 
necessarily forms a caustic. A corollary is that, in cold dark matter cosmology, galactic halos have 
inner caustics in addition to the more obvious outer caustics. The outer caustics are fold catastro- 
phes located on topological spheres surrounding the galaxy. To obtain the catastrophe structure 
of the inner caustics, we simulate the infall of cold collisionless particles in a fixed gravitational 
potential. The structure of inner caustics depends on the angular momentum distribution of the 
infalling particles. We confirm a previous result that the inner caustic is a "tricusp ring" when 
the initial velocity field is dominated by net overall rotation. A tricusp ring is a closed tube whose 
cross section is a section of an elliptic umbilic catastrophe. However, tidal torque theory predicts 
that the initial velocity field is irrotational. For irrotational initial velocity fields, we find the inner 
caustic to have a tent-like structure which we describe in detail in terms of the known catastrophes. 
We also show how the tent caustic transforms into a tricusp ring when a rotational component is 
added to the initial velocity field. 



PACS numbers: 98.80 Cq 
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I. INTRODUCTION 



It is generally believed at present jjj that approximately 23% of the energy density of 
the universe is in the form of "cold dark matter" (CDM). The CDM particles must be non- 
baryonic, collisionless, and of small primordial velocity dispersion. Particle physicists have 
put forth several candidates with the required properties. Two among these, the axion and 
the neutralino, have the distinction of having been originally postulated for purely particle 
physics reasons. The axion solves the "strong CP problem", whereas the neutralino is a 
prediction of supersymmetric extensions of the Standard Model. The primordial velocity 
dispersions of these two candidates are very small, of order 3 • 10~ 17 c for axions and 10~ 12 c 
for neutralinos. 

A central problem in dark matter studies is the question of how CDM is distributed in 
the halos of galaxies, and in particular in the halos of spiral galaxies such as our own Milky 
Way. Indeed knowledge of this distribution is essential for understanding galactic dynamics 
and for predicting signals in direct and indirect searches for dark matter on Earth. 

Three main approaches towards determining the dark matter distribution in galactic 
halos have been put forth. The first assumes that galactic halos are isothermal |2|. This 
assumption is highly predictive and, in fact, some of the predictions of the isothermal model 
agree with observation, to wit the flatness of rotation curves and the existence of core radii 
|3j]. However, even if galactic halos were thermalized in the past, they do not remain so 
because surrounding dark matter keeps falling onto them |5|. There is no mechanism 
which thermalizes the flows of axions or neutralinos which fall late onto the galaxy. The 
second approach is to carry out iV-body simulations of the halos 0]. This does presumably 
give a correct description if the number N of particles in the simulations is sufficiently 
large. However present simulations have only iV ~ 10 7 and, as a result, suffer from 2-body 
relaxation 0] and from inadequate sampling of phase space |j| . 

The third approach exploits the fact that the CDM particles lie on a 3-dim. hypersurface 
in 6-dim. phase space 0^ H, 0| . This implies that the velocity distribution of dark matter 
particles at every point in physical space is discrete and that there are surfaces in physical 
space, called caustics, where the density of dark matter particles is very large. We have 
recently argued (Hf that these discrete flows and caustics are a robust prediction of cold dark 
matter cosmology. The reader may wish to consult ref. for background information and 
a list of references. 

In the present paper we investigate the catastrophe structure of the caustic formed in a 
flow of cold collisionless dark matter particles falling in and out of a galactic gravitational 
potential well. We call the caustics thus formed the "inner caustics" of the galactic halo. 
Galactic halos also have "outer caustics". An outer caustic occurs near where an outflow 
of dark matter turns around before falling back in (loL ITU. The catastrophe structure of 
outer caustics is simple. They are fold (A 2 ) catastrophes [12] located on topological spheres 
surrounding the galaxy. We will see below that the catastrophe structure of inner caustics 
is relatively more complicated. 

We start off below (Section II) by proving that there always is an inner caustic, i.e. that 
it is impossible for cold collisionless particles to flow from all directions in and out of a region 
without forming a caustic. The proof indicates that the inner caustic occurs near where the 
particles with the most angular momentum are at their distance of closest approach to a 
central point of the region. This suggests that the structure of an inner caustic depends 
mainly on the angular momentum distribution of the infalling particles. 
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In previous work |8[, it was shown by analytical methods that the inner caustic is a 
tricusp ring when the initial velocity field v(f) of the infalling particles is dominated by net 
overall rotation. A tricusp ring is a closed tube whose cross section is a section of an elliptic 
umbilic catastrophe We confirm this result in Section IV. Fig. |2(b)| shows an axially 
symmetric tricusp ring as it appears in our simulations. 

However, as will be shown below (Section III), the leading theor y fo r the origin of the 
angular momentum of galaxies, namely tidal torque theory [3, FLU Il5|. predicts that the 
initial velocity field of dark matter particles is a pure gradient. A pure gradient field is 
of course irrotational (V x v — 0). It can nonetheless carry net angular momentum L = 
J d 3 r p(r) r x v(f) because the density p is inhomogeneous. 

Motivated by tidal torque theory, we want to determine the structure of the inner caustics 
of galactic halos when the initial velocity field is irrotational and, more generally, when the 
initial velocity field is not dominated by net overall rotation. We address this issue by 
simulating the infall of cold collisionless particles in a fixed gravitational potential well. The 
inner caustic is revealed by finding the locus of points where the Jacobian of the map between 
the initial and final positions of the particles vanishes, i.e. where the map is singular. 

We consider all initial velocity fields which are linear in position r, i.e. fields of the type 
v = Mr where M is a 3 x 3 real matrix. We set Tr(M) = and more generally ignore 
the radial component v r = f ■ v because it does not contribute to angular momentum and 
hence does not have much influence upon the inner caustics. Tidal torque theory predicts 
M T = M, but we consider the more general case in which M has both symmetric and 
anti-symmetric parts. 

The outline of our paper is as follows. In Section |H1 we introduce the formalism of zero 
velocity dispersion flows and caustics, and give an existence proof of inner caustics. In 
Section IIII1 we discuss our numerical techniques and the initial conditions expected from 
tidal torque theory. In Section IV, we simulate flows which are dominated by net overall 
rotation, confirm that the inner caustic is a tricusp ring in that case, and demonstrate the 
stability of the tricusp ring under small perturbations. In Section|Swe simulate flows which 
are not dominated by overall rotation and derive the structure of inner caustics for those 
cases. In Section VI, we summarize our conclusions. 

In this paper, unless stated otherwise, the words "circle" and "sphere" will be used in 
their topological, rather than geometrical, sense. So an ellipse will be called a circle, etc. 



II. EXISTENCE PROOF OF INNER CAUSTICS 

Consider a flow of cold collisionless particles falling from all directions in and out of 
the gravitational potential well of a galaxy. That the particles are collisionless means that 
they move under the influence of purely gravitational forces. We assume in this paper that 
Newtonian gravity applies, and so the particles obey the equation of motion 

% = - V4(r, ( ) (1) 

where $(r, t) is the gravitational potential of the galaxy. Note however that the existence 
proof of inner caustics given below would hold equally well in General Relativity. That the 
particles are cold means that they have negligible velocity dispersion. We set the velocity 
dispersion equal to zero. However, the presence of a small velocity dispersion does not affect 
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the existence of caustics, and does not change their structure. It only cuts off the divergence 
of the dark matter density at the location of the caustics. 

As was mentioned in the Introduction, CDM particles lie on a 3-dim. hypersurface in 6- 
dim. phase space. Because the number of particles is huge - of order 10 84 axions and/or 10 68 
neutralinos per galactic halo - the particles can be labeled by three continuous parameters 
a = (ai, 0:2,0:3) which are chosen arbitrarily. For example, one may choose a to be the 
position q of the particle at some early initial time t in , before shell crossings have occurred. 
Other parametrizations may be more convenient however, depending on the problem at 
hand. 

Let x(a; t) be the position of the particle labeled a, at time t. After shell crossings have 
occurred, there will in general be particles with different values of a at the same location r. 
Let dj(r,t), with j = l...n(r,t), be the solutions of r = x(a;t). Thus n(r,t) is the number 
of distinct flows at location r at time t. n is always a positive odd integer because n = 1 to 
start with (before shell crossings have occurred) and the number of solutions of r = x(a;t) 
can only change by two at a time. Let ^r(o) be the number density of particles in the 
chosen parameter space. The mass density of particles in physical space is then pj 



The magnitude of D is the Jacobian of the map a — > x. Caustics occur where D — 0, 
i.e. where the map is singular. Note that, although D is not reparametrization invariant, 
the density d and the caustic condition D = are reparametrization invariant. Caustics lie 
generically on 2-dim. surfaces because, in general, the condition D = defines a surface in 
physical space. Only in special cases does the condition D = define an isolated line or 
point. Hence isolated line caustics and isolated point caustics are degenerate cases; they are 
unstable towards becoming caustic surfaces. Finally note that the map a — > x is singular 
where the number of flows n changes. So caustics lie generically at the boundaries between 
regions which have different numbers of flows. On one side of a caustic surface are two more 
flows than on the other. 

We now show that a continous flow of CDM particles falling in and out of a gravitational 
potential well cannot be free of caustics. Let us define a geometrical sphere of radius R 
surrounding the potential well. The precise value of R and the precise location of the 
sphere's center do not matter. Let us label the particles by a = (8, <p, r) where r is the time 
when the particle crosses the sphere on its way into the well, and 9 and are the polar 
coordinates of its crossing point on the sphere. x(9, 0, r; t) is the position of the particle as 
a function of time t. We will show that, at any t, 



vanishes at at least one point inside the sphere. The flow is considered at an arbitrary fixed 
time t. We will suppress the label t henceforth. 




(2) 



where m is the particle mass and 




(3) 
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For each (#, 0), the initial crossing time parameter r has some range: T ou t(0, 0) < r < 
7i n (0, 0) where T m (r out ) is the initial crossing time of particles presently crossing the sphere 
on the way in (out). [The careful reader may notice that Ti n (9,(p) = t, but this does not 
play a special role in what's to follow.] Let us choose the origin x = at the center of the 

sphere and consider the function of three variables r(9,(fi, r) = ^Jx(9, 0, r) ■ x(9, 0, r). We 

have ^\e,4,,Tou t (6,<t>) < and Wf\e,<j>,rU0,<t>) > °- Hence for all (0, 0) there exists r (9, 0) 
such that 

r(9, 0, tq(9, 0)) = min r(9, 

where the minimum is over r G [r ont (9, 0), r in (9, 
distance to the origin among all particles labeled 

for some (9, 0) . We will return later to the opposite case, where r 



r) = r min {U, <p) 

)} for fixed (9,(f>). r min (9, 
Kd>). Let us first assume that r 



(5) 



is the closest 

min(0, 0) 7^ 

(9, 0) = for all 



We have 



dr 



for all 



such that TV 



x dx 
6,(/>,T (6,<p) r dr 
,) ^ 0. Now, r min ( 



S-2 = {(9,(j))}. Let (6*0,00) be such that 



= (6) 
, 0) has a maximum value over the sphere 



max r r , 



(7) 



We have then 



dr 
89 



"0 



r 



dx 
~89 



dr 



"0 



"0 



r 



9x 



with do = (9q, 0o, ^(^o, 0o))- Note that = x(do) ^ 0. Eqs. 



I (8) 
and (JEI) imply that 



|f |<3 > ||lao anc ^ ffl«o are an perpendicular to x . Hence those three vectors are linearly 
dependent, and therefore D(a Q ) = 0. This implies that x$ is the location of a caustic. Note 
that xq depends on the choice of origin. If we move the origin about, xq will move too. 
Thus the inner caustic is in general spatially extended. This is as expected since caustics lie 
generically on surfaces. 

Next we consider what happens in the special case where r m i n (#,0) = for all (#,0). 
Then x(9, 0, t (9, 0)) = for all {9, 0). Therefore, for r near t (9, 0), 



x(9,<P,t) = v(9, 



ro(^,0)) + O((r-r o ^,0)) 2 ) 



(9) 



where v( 



rename 
t. Hence 



dx 
8t 



3' A' 



. We may reparametrize 9' = 9, 0' = 0, r' = r — tq(9, 0), and 

4>,t (o,4>) 

t') — > (9, 0, r). In this parametrization, x(9, 0, r) = u{9, 0)r + O(r 2 ) for small 



(10) 



Since D = at r = 0, the origin is the location of a caustic. This completes the proof. Note 
that the case xo = is special since the whole inner caustic has collapsed to a point. 

To conclude this section, let us remark that in the situation discussed here, where CDM 
particles fall in and out of a gravitational potential well, the number of participating flows 
inside the sphere with radius R is everywhere an even integer. Indeed just inside the r = R 
surface there are two flows, one going in and one going out, and (as remarked earlier) the 
number of flows can only change by two at a time. The reader may wonder how this fits 
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with the statement that the number of flows is everywhere odd. The resolution of this 
little puzzle is that, in actual realizations such as a galactic halo, there is always an odd 
number of additional flows present which are not participating in the in and out flow under 
consideration. 



III. INITIAL CONDITIONS AND NUMERICAL TECHNIQUES 

In this section we describe how the simulations are done. A central issue is the initial 
conditions we give to the particles. In subsection A we show that in the tidal torque theory 
for the origin of the angular momentum of spiral galaxies the flow of CDM particles is 
irrotational to all orders in perturbation theory, i.e. to all orders in an expansion in powers 
of the density perturbations. In subsection B, we present the initial conditions used in the 
simulations. They are a generalization of the prediction of tidal torque theory to allow for 
rotational flow as well as irrotational flow. In subsection C, we describe the steps carried 
out in doing the simulations. 



A. Irrotational flow in tidal torque theory 

Net rotation is a striking property of isolated spiral galaxies. Yet it is not clear at present 
that the origin of this net rotation is well understood [l6|. The leading hypothesis is that net 
rotation of spiral galaxies is the result of torque applied by the tidal gravitational forces of 
neighboring density perturbations in the very early stages of structure formation [l3l . 13 . ljf . 



This hypothesis is called "tidal torque theory". It is assumed in this context that general 
relativisitic effects are unimportant, and that Newtonian gravity applies. 

In CDM cosmology, density perturbations enter the non-linear regime when shell crossings 
occur and caustics form. Before that time, or wherever shell crossings have not occurred 
yet, there is a single flow at every physical location f. Let us call v(f, t) the velocity of that 
primordial flow. In the formalism of Section II, v(f, t) is obtained by eliminating a from 
f= x(a,t) and v — |f (a, t). For collisionless particles 

dv dv —* — * 

= -Ql^t) + (v(f,t) ■ V)v(r,t) = -V$(f,t) (11) 

where $(r, t) is the gravitational potential. Eq. (fTT]) neglects relativistic effects such as the 
gravitomagnetic force of General Relativity. Eq. (fTTjl implies 

d - 

_(V x v) = -V x [(v-V)v\ 

= -itijk[(djVi)(diVk) +vididjV k ] . (12) 

If v = V%, both terms on the RHS of Eq. ()12j) vanish. Hence a pure gradient initial velocity 
field remains pure gradient at all times. 

In zeroth order of perturbation theory, the flow is given by Hubble's law: v(f,t) = H(t)f 
where H(t) — - and a(t) is the scale factor. That velocity field is certainly a pure gradient. 
In first order, the particle trajectories are given by 

x(q,t) = a{t)[q-b{t)V q <S>{Q)\ (13) 
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where q is the particle's position at a very early time, and <&(q) is the gravitational potential. 
[The latter is time independent in first order perturbation theory when expressed in terms 
of the co-moving coordinate q, i.e. $(r = a(t)q,t) = $(<?*).] In the matter dominated era, 
both a(t) and b(t) are proportional to ts. As expected, the velocity field implied by Eq. (JT3Jl 
is irrotational. Our remark implies that the velocity field remains irrotational to all orders 
of perturbation theory. 

B. The linear initial velocity field approximation 
Eq. (JT3j) implies 

db 

v(r,t) = H(t)f- a(t)j t (t) V q m\ q -=^r • (14) 

Following ref. let us choose q = at a minimum of $ and expand $ in Taylor series 
up to second order in the q's. The velocity field is then a linear function of position: 

v{f) = M r (15) 

where M is a symmetric (M T = M) 3x3 matrix. When a shell of particles surrounding 
a galaxy approaches turnaround, its initial Hubble expansion has been canceled by the 
attractive gravitational force of the galaxy. For such a shell the trace of M is approximately 
zero. 

In Sections IV and V we numerically integrate the equations of motion of particles falling 
in and out of the galactic gravitational potential well. The particles are given initial condi- 
tions on a geometrical turnaround sphere as follows. The particle labeled (9, <f>, r) starts at 
time r at location 

n a (9,4>) = Rn(9,4>) (16) 

where n(9, (f>) is the unit vector in the direction defined by polar angle 9 and azimuth </>. R 
is the radius of the turnaround sphere. The particles are assigned velocities of the form of 
Eq. (USD with Tr M = 0. We write M = S + A with S T = S and A T = -A. Tidal torque 
theory predicts A = 0. However, as tidal torque theory may not be the final word on the 
origin of the angular momentum of galaxies, we will study the inner caustics for A ^ as 
well as A = 0. 

We believe the above set of initial conditions is sufficiently general for our purpose of 
studying the structure of inner caustics. Indeed, the structure of inner caustics is determined 
by the distribution of distances of closest approach of the particles falling in. This in turn is 
determined by the distribution of angular momenta on the turnaround sphere, or equivalently 
by the tangential components of the initial velocity field on the turnaround sphere. Caustics 
are stable under deformations. So to find the structure of inner caustics it is not necessary 
to use the most general initial conditions, but only initial conditions which are sufficiently 
representative of those that occur in reality. We expect that to include higher order terms 
in the Taylor expansion of the initial velocity field as a function of position will only deform, 
without changing their essential structure, the inner caustics found when assuming the initial 
velocity field is of the form of Eq. ()15|) . We verify this expectation explicitly in the case of 
the tricusp caustic ring when we study its stability in Section IV. 
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We choose coordinate axes such that S is diagonal: 



M 



R 



1 



9i -C3 c 2 

C 3 #2 -Ci 

-c 2 ci -g x - g 2 



) 



(17) 



gi and #2 parametrize the symmetric part of M which yields the gradient part of v, whereas 
Cj, c 2 and C3 parametrize the anti-symmetric part of M which yields the curl part of v and 
describes a rigid rotation of angular velocity uj = -|. In terms of the five parameters, the 
components of the initial velocity field tangent to the turnaround sphere are 



= (#2 — gi) sin 9 sin cos — cos 6{c\ cos + C2 sin 0) + C3 sin 9 
v e {9,4>) = v-6 

= sin9 cos 9 [g%(l + cos 2 0) + g 2 (l + sin 2 0)] — ci sin + c 2 cos . (18) 

The radial component v ■ f of the initial velocity field does not contribute to the angular 
momentum, and is set equal to zero. We verify in Section IV that the inclusion of radial 
initial velocities on the turnaround sphere has very little effect on the position of tricusp 
caustic rings, and no effect on their structure. 

To conclude this subsection, we discuss the symmetry properties of our initial velocity 
field. Almost always we take the gravitational potential to be spherically symmetric, in 
which case the symmetry properties of the initial velocity field are those of the subsequent 
evolution as well. In the irrotational case (ci = C2 = C3 = 0), the initial velocity field is 
reflection symmetric about the x — 0, y — and z = planes. Moreover it is axially 
symmetric when two of the three eigenvalues (gi, g%, and g 3 = —g\ — g 2 ) are equal. Most 
often we will chose the axes such that g\ < g 2 < g%. The parameter space is then g\ < 
and g\ < g% < ~\g\- When g\ = g 2 , the initial velocity field is axially symmetric about the 
z axis. When g 2 = —hgi, it is symmetric about the x axis. 

In the case of pure rotation (gi = g 2 = 0), we may choose axes such that c — ci. The 
initial velocity field is always axially symmetric in this case. When g±, g 2 , c%, c 2 and C3 are 
all different from zero, the initial velocity field has no symmetry. When axial symmetry 
about the z axis is imposed, c\ = c 2 = and g\ = g 2 . 

C. How the simulations are done 

We simulate a single flow of zero velocity dispersion falling in and out of a time inde- 
pendent gravitational potential $(r), which is specified below. The initial conditions are 
Eqs. ()16I18|) plus v r = v ■ f = 0. We solve the equation of motion, Eq. (JTJ), numerically to 
obtain for all (9, 0) the trajectory x(9, 0, r; t) of the particle which started at position (9, 0) 
on the turnaround sphere, at time r. Since neither the potential $ nor the initial conditions 
are time-varying, the simulated flows are stationary, i.e. x(9,(f>,r;t) = x(9,(f),t — r). The 
simulation of non-stationary flows would be straightforward but considerably more mem- 
ory intensive and time consuming, without being more revealing of the structure of inner 
caustics. The only change with respect to simulations of stationary flows would be that the 
caustics deform in a time dependent way. This would not teach us anything new about the 
structure of inner caustics. 



0) 
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From x(9,(p,t' = t — r) we calculate D = det f/f^~f) • We then plot the points where 
D — 0. The set of these points is the inner caustic surface. 

Unless stated otherwise, $ is the gravitational potential produced by the matter density 
profile: 



v 2 



* r » = «J(#T^) • (19) 

which implies an asymptotically flat rotation curve with rotation velocity v TO t- a is the core 
radius. We will refer to the density profile of Eq. (|19|) as the "isothermal" profile although 
it is only an appoximation to an exact isothermal profile. The equation of motion is then 



at 2 r L r \a 



(20) 



We use everywhere R, the radius of the turnaround sphere, as the unit of distance and v rot 
as the unit of velocity. So our unit of time is — , which is of order the infall time. The core 

J Vrot ' 

radius a is always set equal to 0.0285. 

The five parameters gi, g 2 , Ci, c 2 and c 3 , when expressed in units of v Tot , are related to 
and of order the dimensionless angular momentum j defined in ^| . It was estimated in that 
paper that the average of j over the turnaround sphere is approximately 0.2 in the case of 
the Milky Way. This sets the overall scale for the values of gi ... C3 we are interested in and 
which are used in our simulations. Note that it is the relative values of these five parameters 
that determines the structure of inner caustics. The overall scale of the parameters merely 
determines the overall size of the inner caustics relative to R. 

Since we simulate the infall of a single cold flow in a fixed potential, the particle resolution 
is not a critical issue. We chose a resolution of one particle per degree interval in 9 and <p, 
and a time step of 10 -4 . 



IV. THE TRICUSP RING 

It was found in ref. [8| that the inner caustic is a "tricusp ring" when the initial velocity 
distribution is dominated by a rotational component. A "tricusp ring" is a closed tube whose 
cross section has the tricusp shape characteristic of the elliptic umbilic {D-4) catastrophe. 
In this section, we perform simulations of flows with initial velocity fields dominated by a 
rotational component, show that the inner caustic is indeed a tricusp ring and that this 
structure is stable under perturbations. 



A. Axially symmetric case 

Consider the simple case where the initial velocity field is a linear function of the coor- 
dinates and is purely rotational, i.e. when Eq. (JT5)l holds with an anti-symmetric matrix 
M. The initial velocity field is then axially symmetric. Let us choose coordinates such that 
v = C3 sin 6 (f). Fig. Q shows the resulting infall of a single shell in xz cross section, at 
successive times, for c 3 = —0.1. 

As the shell falls in, deviations from spherical symmetry appear due to the presence of 
angular momentum. The particles at the poles have zero angular momentum and fall in faster 
than the particles near the equator [l(b)|. They are the first to cross the xy plane [l(c)|. 
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The shell completes the process of turning itself inside out in Fig. |l(e)[ forming a crease. 
Finally the shell increases in size again and regains an approximately spherical shape |l(f)| . 

The inner caustic occurs at and near the location of the crease in Fig.CJ Fig. |2(a71 shows 
the flow near the crease in pz cross section where p = y/x 2 + y' 2 . The figure shows that the 
pz plane is divided into two regions, one with two flows at each point and the other with 
four flows at each point. The boundary that separates these two regions is the caustic. The 
dark matter density is infinite there in the limit of zero velocity dispersion. Fig. |2(b)| shows 
the caustic in three dimensions. It is indeed the "tricusp ring" described in ref. 8]. 

B. Perturbing the initial velocity field 

1. Breaking axial symmetry 

We introduce departures from axial symmetry by adding gradient terms, proportional 
to gi and g 2 - Fig- El shows the inner caustic for the initial velocity field of Eqs. (fTSj) with 
c 3 = —0.1, gi = —0.033, g 2 = 0.0267. It is again a tricusp ring but the cross section is now 
0-dependent. The tricusp shrinks to a point four times along the ri ng. In the neighborhood 
of each such point, the catastrophe is the full elliptic umbilic(D L 4 ) [12]. 

2. A random perturbation 

We added randomly chosen perturbations to the previously discussed axially symmetric 
initial velocity field. The inner caustic shown in Fig. E]is for the initial velocity field 

v = —0.1 sin8 cf) 

- 0.01 sin 6 [sin(20) + 0.5 sin(30) + 0.25 sin(40)] (21) 

Although the caustic is deformed from what it was in Fig. El its structure remains a ring 
whose cross section is everywhere a tricusp. 

3. The effect of radial velocities 

We added radial velocities to the previously discussed axially symmetric initial velocity 
field. We find that the radial velocity components result in only relatively small changes in 
the dimensions of the tricusp ring. For the initial velocity field 

v = c 3 sin6 ((f) + f) (22) 

with C3 = —0.1, the tricusp ring radius was decreased by 0.28%, compared to what it was 
for the original initial velocity field (v = C3 sin 9 <fi) , and the transverse dimensions of the 
tricusp ring were reduced by 11% and 14% in the directions perpendicular and parallel to 
the plane of the ring. For C3 = —0.3, the radius was reduced by 1.7% and the transverse 
dimensions by 69% and 73% respectively. 

Let us explain why radial velocities on the turnaround sphere have only a small effect 
on the inner caustics. The inner caustics are determined by the distribution of distances 
r min of closest approach to the galactic center of the infalling particles. The distance of 
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closest approach is determined by angular momentum conservation: £ = r min v mSLX} where £ 
is specific angular momentum and v max is the speed at the moment of closest approach. The 
latter is determined by energy conservation 

^Lx = \H + v 2 e + vl) + - $(r min ) . (23) 

The main contribution to v max is from the gravitational potential energy released while the 
particle falls in. The initial velocity components provide only corrections to f max which are 
second order in v^, vg and v r . Since £ does not depend on v r at all, radial velocities produce 
only second order corrections to the distances of closest approach. 



C. Modifying the gravitational potential 

In this subsection we verify that, when the initial velocity field is dominated by a rota- 
tional component, the inner caustic is a tricusp ring independently of the choice of gravita- 
tional potential. 



1. The NFW density profile 



We carried out simulations of the infall of collisionless particles in the gravitational po- 
tential produced by the density profile of Navarro, Frenk and White ET 



p(r) 



Ps 



(24) 



The scale length r s was chosen to be 25 kpc. p s was determined by requiring that the 
rotational velocity at galactocentric distance r = 8.5 kpc is 220 km/s. The acceleration of 
a particle orbiting in the potential produced by the NFW density profile is then: 



a(r 



(220 km/s) 



x 



r x z 



ln(l 



x 



x 
1+x 



ln(l + z ) - 



l+x 



(25) 



where x = r/r s and x Q = r Q /r s . 

Fig. El shows the result of two simulations plotted on the same figure. The outer caustic 
ring is obtained using the density profile of Eq. (|24jh while the inner caustic ring is obtained 
using the density profile of Eq. (|T§|) with t> rot = 220 km/s and a = 4.96 kpc. In both cases, 
the turnaround radius R = 174 kpc and the initial velocity field v = 0.2 sin# (p. As always, 
the coordinates x, y, z are in units of R. The inner caustic is a tricusp ring in each case, but 
with different dimensions. The ring caustic produced by the NFW profile has larger radius 
than that produced by the isothermal profile because the NFW gravitational potential is 
shallower than the isothermal one at the location of the caustic (r caust i C — 16 kpc). Since 
£ = r m i n f max is the same, r min is larger in the NFW case because t> max is smaller. 
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2. Breaking spherical symmetry 



We also simulated the infall of collisionless particles in a non spherically symmetric grav- 
itational potential. For the latter, we chose the triaxial form: 



$fr) = - 



R 



\ 



i) + (i) + te) ) 



(26) 



where a 1; a 2 and 03 are dimensionless numbers. Fig. |U] shows the inner caustic for the case 
where a% = 0.95, a 2 = 1.0 and a 3 = 1.05, and the initial velocity field v = 0.2 sin^ cj). It is 
again a tricusp ring. Its axial symmetry is lost due to the absence of axial symmetry in the 
potential. The tricusp ring still has reflection symmetry about the xy, yz and xz planes. 
As in Fig. |H] the tricusp shrinks to a point four times along the ring. 



V. GENERAL STRUCTURE OF INNER CAUSTICS 

In this section, we describe the structure of inner caustics when the initial velocity field is 
not dominated by a rotational component. In subsection A, we discuss the axially symmetric 
case, whereas the non axially symmetric case is discussed in subsection B. In each subsection, 
we simulate first irrotational (i.e. pure gradient) initial velocity fields. As was mentioned in 
Section III, irrotational initial velocity fields are predicted by tidal torque theory. We will 
find the inner caustics produced by irrotational velocity fields to have a definite structure 
which we refer to as the "tent caustic" . After describing the tent caustic, we add a rotational 
component to the initial velocity field and see how the tent caustic deforms into a tricusp 
ring. 



A. The axially symmetric case 

The initial velocity field of Eqs. (|18|) is symmetric about the z axis when c\ = c 2 = and 
gi = g 2 . Then 

3 

v = ~gi sin(2#) 9 + c 3 sin 9 <j> . (27) 

We first simulate the flow and obtain the inner caustic in the irrotational case, C3 = 0. Next 
we see how the flow and inner caustic are modified when c 3 ^ 0. 



1. Infall of a cold collisionless shell 

Figs. [7| and |H1 show the infall of a cold collisionless shell whose initial velocity field is 
given by Eq. (|27j) with c 3 = and g\ = —0.0333. Since c 3 = 0, each particle stays in the 
plane containing the z axis and its initial position on the turnaround sphere. Figs. and 
IH1 show the particles in the y = plane. The angular momentum vanishes at 9 = and 
9 = 7r/2 where 9 is the polar coordinate of the particle at its initial position. Hence, the 
particles labeled 9 = and 9 = 71/ 2 follow radial orbits. The angular momentum increases 
in magnitude from 9 = tt/2, reaches a maximum at 9 = ir/A and returns to zero at 9 = 0. 
The sign of the angular momentum does not change during this interval. 
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The shell starts out as shown in Fig. |7(a)] As the shell falls in, the particles at 9 7^ 0, ir/2 
move towards the poles. These particles feel an angular momentum barrier and fall in more 
slowly than the particles at 9 = 0. This results in the formation of a loop in Fig. |7(c)| 
The formation of the loop implies a cusp caustic on the z axis. The particles labeled 9 = 
and 9 = 7T have crossed the z = plane and the particles labeled 9 = tt/2 have crossed the 
x = plane in Fig. |7(g)| The shell then takes the form shown in Fig. |7(h)| The further 
evolution is shown in Fig. |S| The loop that is present near the z = plane in Fig. 8(b) 



disappears through the sequence of Figs. |8(c)| - |8(g)[ The disappearance of the loop implies 
the existence of a cusp caustic in the z = plane as well. In Fig. |8(hj] the shell has regained 
an approximately spherical form and is expanding to its original size. 

For larger values of |gi| the early evolution is qualitatively the same as in Fig. [71 but the 
late evolution is qualitatively different from Fig. |HJ Fig. El shows the late evolution of a cold 
collisionless shell with the initial velocity field of Eq. (J27j) with C3 = and g\ = —0.0667. 
The loop which is present near the z 
sequence 



9(c) - 9 g 



plane in Fig. |9(b) 
than was the case in Fig. [HI In Fig. 



the z = plane before the sphere turns itself inside out 
structure, and a more complicated caustic, than for the g\ 



disappears by a more complicated 
the particles near 9 = n/2 cross 
This crossover produces additional 
—0.0333 case. The critical 



value of |<7i|, below which the qualitative evolution is that of Fig- [HI and above which that 
of Fig. El is g u ~ 0.05. 



2. Caustic Structure 



The inner caustic is a surface of revolution whose cross section is shown in Fig. |10(a)] for 
the case (03,(71) = (0, —0.0333) and in Fig. |11 (ai)] for (03,(71) = (0, —0.0667). For the sake of 
brevity, we call this structure a "tent caustic". On the z axis, there is a caustic line which 
we call the "tent pole". The remainder of the caustic is called the "tent roof. As was 
mentioned in Section II, caustic lines are not generic. The tent pole is a line in Figs. |10(a)| 
and |ll(a)| only because the initial velocity field is axially symmetric and irrotational. We 
will see below that when axial symmetry is broken or when a rotational component is added, 
the tent pole becomes a caustic tube of a specific sort. 

For \gx\ < gu there is a cusp in the tent roof where it meets the z = plane and two 
cusps where the tent roof meets the tent pole, one at the top and one at the bottom. For 
gi < —gu the cusp in the z = plane is replaced by a butterfly catastrophe The 
butterfly has three cusps and three points of self-intersection. The cusp in the z = plane 
transforms into a butterfly by increasing \gi\. The latter parameter may therefore be called 
the "butterfly factor" Q. 

If gi is chosen positive instead of negative, the behavior at the poles and the equator 
is reversed [see Eq. (|27|)] and we have cusps on the x axis and either cusp or butterfly 
catastrophes on the z axis, depending on the magnitude of g±. 

Fig. |1 0(b)! shows the dark matter flows in the vicinity of the caustic for \gi\ < gu. There 
are four flows everywhere inside the caustic tent and two flows everywhere outside. Fig. |1 1(b) 
shows the dark matter flows in the vicinity of the butterfly caustic, for g\ < —gu- Fig. 112 
shows the number of flows in each region of the butterfly caustic. 
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3. Adding a rotational component 



Here we show, in the axially symmetric case, the effect of adding a rotational component 
to the initial velocity field. On the basis of the discussion in Section IV, we expect the tent 
caustic to transform into a tricusp ring. Fig. IT3*1 shows the transformation. We start with an 
irrotational velocity field (03 = 0) in 13(a)| and increase C3 until the rotational component 
dominates the velocity field, in |13(d) 



The tent pole on the z axis, which is a line caustic 
in the irrotational case, changes to a tube of circular cross section. The radius of this tube 
increases until the tent caustic becomes indistinguishable from a tricusp ring. 



B. Non axially symmetric case 

We start off by discussing the flows and caustics resulting from irrotational initial velocity 
fields. With c = 0, Eqs. ([T%]) become 

v = £ sin6 sin(20) + sin(2#) (f sin 2 + g) 9 , (28) 

where £ = \{§2 — Qi) and g = g\ + \gi- £ is a measure of z axial symmetry breaking in 
the irrotational case. We first let £ <C g. Next, we explore all of (<?i,<72) parameter space. 
Finally, we add a rotational component by letting C3 7^ 0. 



1. Irrotational, non axially symmetric perturbations 

We saw in subsection V-A that there is a caustic line on the z axis (the tent pole) when the 
initial velocity field is irrotational and axially symmetric. Fig. |14(a)| shows the trajectories 
of the particles in the z = plane for such a case. The orbits are radial. Indeed all particles 
have zero angular momentum with respect to the z axis when the initial velocity field is 
irrotational and axially symmetric. Because all trajectories intersect the z axis, there is a 
pile up of particles on that axis and hence a caustic line. 

Fig. 14(b) I shows the trajectories of the particles in the z = plane for the initial velocity 



field of Eq. ()28|) with £ = 0.01 and g = —0.05. The particles do have angular momentum 
with respect to the z axis now. The caustic line on the z axis spreads onto a tube whose 
cross section is the diamond shaped envelope of particle trajectories shown in Fig. |14(b)[ 
That envelope has four cusps. The flows and caustic have reflection symmetry about the 
xy, xz and yz planes because the initial velocity field has those symmetries. Fig. |14(b)| 
shows four flows inside the diamond-shaped caustic and two flows outside. The infall of 
dark matter particles from regions above and below the z = plane will add two more flows 
at each point, which are not shown in Fig. El for clarity. 

Fig. El shows the inner caustic in 3 dimensions for the initial velocity field of Eq. (J28|) 
with £ = 0.005 and g = —0.05. The tent pole has spread onto a tube with diamond shaped 
cross section, as in Fig. [j~4(bj| Fig. [15(b)] shows a succession of constant z sections. The 
(topological) circles, which are sections of the tent roof, and the diamond structure can be 
seen clearly. Near the z = plane, there are six flows inside the diamond, four flows in the 
other regions inside the caustic tent, and two outside the tent. Fig. El shows y = and 
x = sections of the tent caustic. 
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2. A hyperbolic umbilic catastrophe 



Let us look more closely at the two regions (top and bottom) in Fig. EH where the tent 
pole reaches and traverses the tent roof. Figs. |17(a)| - |17(d)| show z = constant sections of 
the inner caustic in such a region. As \z\ is increased, the two cusps of the pole which are in 
the x = plane simply pierce through the roof, whereas the parts of the pole near the y = 
plane traverse the roof by forming with the latter two hyperbolic umbilic (D +4 ) catastrophes 
|l2| , one on the positive x side and one on the negative x side. The sequence through which 
this happens is shown in greater detail in Figs. |17(e)| - |17(h)| for the hyperbolic umbilic on 
the positive x side. The arc (section of the tent roof) and the cusp (section of the tent pole) 
approach each other until they overlap [17(g)| , forming a corner. After they have crossed, 
the roof section is cuspy whereas the pole section is smooth. The cusp is transferred from 
the tent pole to the tent roof as the two surfaces pass through one another. This behaviour 
is characteristic of the hyperbolic umbilic catastrophe. There are four hyperbolic umbilics 
embedded in the caustic tent structure, two (x > and x < 0) at the top (z > 0) and 
two at the bottom (z < 0). The hyperbolic umbilic at z > and x < is shown in three 
dimensions in Fig. |17(i)| 

Let's mention that the particles forming the pole and roof where they intersect in the 
x = plane originate from different patches of the initial turnarround sphere whereas the 
particles forming the pole and roof near a hyperbolic umbilic originate from the same patch 
of the initial turnaround sphere. 



3. The (51,(72) landscape 

Here we describe the inner caustic in the irrotational case for g\ and g 2 far from those 
values where the flow is axially symmetric. Recall that the flow is symmetric about the z 
axis when g 1 = g 2 (#3 = — 2t7 X ) , about the y axis when g 2 = —2g 1 (g 3 = gi), and about the 
x axis when g 2 = —\g\ (#3 = g2)- in terms of £ and g, these conditions for axial symmetry 
are £ = 0, g = 0, and £ = — g, respectively. 

The first, second and third column of Fig. El show respecively the z = 0, y = and 
x = sections of the inner caustic produced by the initial velocity field of Eq. (|2~%)) for 
various values of (pi,^)- The ratio g 2 /g\ decreases uniformly from 1 (top row) to -1/2 
(bottom row). Note that the third column describes a sequence which is that of the first 
column in reverse, and that the first half of the sequence in the second column is the reverse 
of the sequence in its second half, with x and z axes interchanged. 

In the first row, the caustic is axially symmetric about the z axis. It is as described earlier 
in Fig. In the second row, the axial symmetry is broken and the tent pole has acquired 
a diamond shaped cross section. The caustic is now as described in Figs. [IS] - H3 In the 
third row, the cusps on the tent pole that lie in the x = plane have pierced the tent roof 
all the way from top to bottom, and the hyperbolic umbilics in the y = plane have moved 
towards the z = plane. In the fourth row, the hyperbolic umbilics have almost reached 
the z = plane. What was the tent roof in the first row is now stretched along the x axis, 
and becomes the tent pole in the fifth row when the inner caustic is axially symmetric about 
the x axis. In the process described by Fig. EE i n which the inner caustic metamorphoses 
from a tent symmetric about the z axis to a tent symmetric about the x axis, the tent roof 
smoothy deforms into the tent pole and vice-versa. 
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The plots of Fig. EH are reminiscent of caustics seen in gravitational lensing theory 21 
221 . 23^ . and ship stability analysis 2^, 25^ . 



4- Adding a rotational component 

Here we add a rotational component (03 7^ 0) to the initial velocity field of Eq. ()28|) to see 
the metamorphosis of the inner caustic from tent to tricusp ring. Fig. shows the z = 
cross sections of the inner caustic during such a transition. In Fig. |19(a)[ the i nitial velocity 
field is irrotational and we see a circle and diamond, as before. In Fig. |19(b)[ the diamond 
is skewed because of the rotation in the z = plane introduced by c 3 7^ 0. Fig. |1 9(c)] shows 
the case C3 = £. As C3 is increased further, the diamond transforms into two swallowtail 
(A4) catastrophes l3] joined back to back 19(d) . There are two flows in the central region, 



six flows in the cusped region of each swallowtail, four in the other regions inside the circle 
and two outside the circle. Finally, the swallowtails pinch off to form the inner circle of the 
z = section of the tricusp ring. Fig. [201 shows the transition in three dimensions. 

VI. CONCLUSIONS 

In Section II, we gave a mathematical proof of the statement that a cold flow of collision- 
less particles from all directions in and out of a region necessarily produces a caustic. We 
call this caustic the "inner caustic" of the in and out flow. The main purpose of our paper 
was to determine the catastrophe structure of the inner caustics formed by cold dark matter 
particles falling in and out of a galactic gravitational potential well. 

The structure of the inner caustic depends for the most part on the angular momentum 
distribution of the particles falling in. It had been shown previously by analytical methods 
[H that the inner caustic is a tricusp ring when the velocity distribution of the infalling 
particles is dominated by net overall rotation. However we show in Section III that the 
leading theory for the angular momentum of galaxies, namely tidal torque theory, predicts 
that the velocity field is irrotational (V X v — 0) to all orders of perturbation theory, i.e. 
to all orders in an expansion in powers of the size of density perturbations. So, tidal torque 
theory states not only that the initial velocity field of the infalling particles is not dominated 
by net overall rotation, it states that the initial velocity field is exactly irrotational. 

So the question was: what is the structure of the inner caustic when the initial velocity 
field is irrotational? Or more generally, what is that structure when the initial velocity field 
is not dominated by net overall rotation? We addressed this issue by simulating the flow of 
cold collisionless particles falling in and out of a fixed gravitational potential. In principle we 
can do this for any initial velocity field. However, we restricted ourselves to initial velocity 
fields of the form v = Mx where x is initial position and M is a 3 x 3 real traceless matrix. 
As far as uncovering the catastrophe structure of inner caustics is concerned, we expect this 
to be a sufficiently broad set of initial velocity fields. Adding higher order terms to the 
expansion of v in powers of x is expected to merely deform the inner caustics obtained when 
keeping linear terms only. M can be written in the form of Eq. (|17|) which depends on five 
parameters: c— (ci,c 2 ,c 3 ), g\ and g 2 . The first three describe net rotation with angular 
velocity c. The last two describe irrotational flow. 

In Section IV, we simulated flows which are dominated by net rotation and confirm that 
the inner caustic is a tricusp ring in that case. We show that the tricusp ring is stable 
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under perturbations both in the initial velocity field and in the gravitational potential. This 
stability is not a surprise since it is known that the structure of non-degenerate catastrophes 
is stable under perturbations j2^|. When the initial velocity field and/or the gravitational 
potential are not axially symmetric, the dimensions of the tricusp vary along the ring. In the 
neighborhood of a point where the tricusp dimensions have shrunk to zero, the catastrophe 
structure is the full elliptic umbilic. 

In Section V, we simulated flows which are not dominated by a rotational component. 
We started off simulating flows which are both irrotational (c = 0) and axially symmetric 
{gi — 92)- For such flows the inner caustic is shown in Fig. ^] for \g\\ < ~ 0.05. We 
call this structure a "tent caustic" . It has a caustic line on the axis of symmetry, which we 
call the "tent pole", connected to a caustic surface which we call the "tent roof. When 
gi < —gx*, the inner caustic is the same as for \g x \ < g^ except the cusp in the equatorial 
plane is replaced by a butterfly caustic. When g 1 > gi*, the cusps on the symmetry axis are 
replaced by butterfly caustics. 

When g 2 is different from but close to g x , axial symmetry is slightly broken. In that case 
the caustic pole spreads onto a tube whose cross section has four cusps forming a diamond 
shape. See Figs. and E3 In a region where the tent pole connects with the tent roof, 
two hyperbolic umbilic catastrophes appear, as described in Fig. [T7| Fig. [THI shows what 
happens when g 2 is very different from g\. The inner caustic metamorphoses from a tent 
which is symmetric about the z axis when g 2 — g\ to a tent symmetric about the x axis 
when g 2 = -\g x . 

Finally, we investigated the smooth transformation of the tent caustic into a tricusp ring 
when a rotational component is added (c^ 0) . In the axially symmetric case (see Fig. ITHJ) 
the tent pole spreads onto a tube of circular cross section. The radius of this tube becomes 
the inner radius of the tricusp ring. So, in the axially symmetric case, we may think of 
the tent caustic as a tricusp ring whose inner radius has shrunk to zero. In the non axially 
symmetric case, the metamorphosis of tent caustic to tricusp ring is shown in Figs. EH and 
|2*U1 The diamond shaped cross section of the tent pole transforms to a (topological) circle 
by passing through an intermediate stage where it consists of two swallowtail catastrophes 
back to back. 
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FIG. 1: Infall of a cold collisionless shell for the axially symmetric initial velocity field v = 
—0.1 sin# (j>. The six frames show the shell in cross section at times t' = 1.206, 1.258, 1.263, 
1.265, 1.267, 1.275. The continuous infall of many such shells produces the caustic shown in 
Figs. 2(a) and 2(b) x and z are in units of the outer turnaround radius R. Note that here, and in 



many subsequent figures, the x and z scales vary from frame to frame. 
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(a) 




(b) 

FIG. 2: (a) Dark matter flows near the inner caustic for the initial velocity field v = —0.1 sin# (p. 
The inner caustic is a ring whose cross section has three cusps, one of which points away from the 
galactic center. The three cusps are clearly visible in the figure. There are four flows everywhere 
inside the caustic surface, and two everywhere outside, (b) The axially symmetric tricusp ring in 
three dimensions. Note that here, and in many subsequent figures, the caustic is stretched in the 
z direction for greater clarity. 
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FIG. 3: Tricusp ring for the initial velocity field of Eq. I)18j) with C3 = —0.1, g\ = —0.033, §2 
0.0267. It has reflection symmetry about the xy, xz and yz planes. 




FIG. 5: Tricusp rings caused by the same initial velocity field but different gravitational potentials. 
The larger ring (green) was obtained using the NFW profile, Eq. (|24|). while the smaller ring (red) 
was obtained using the density profile of Eq. (|19|). 




FIG. 6: Tricusp ring for the initial velocity field v = 0.2 sin# <p and the non spherically symmetric 
gravitational potential of Eq. (|26|) . 
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FIG. 7: Infall of a cold collisionless shell for the axially symmetric initial velocity field v = 
—0.05 sin(2#) 9. The shell is shown in xz cross section at times t' = 1.132, 1.235, 1.241, 1.251, 
1.256, 1.258, 1.260 and 1.261. The further time evolution is shown in Fig. El The continuous infall 



of many such shells produces the caustic shown in Fig. 10(a 
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(g) 00 

FIG. 8: Continuation of the sequence of Fig. The shell is shown at times t' = 1.265, 1.269, 1.271, 
1.272, 1.273, 1.2737, 1.2742 and 1.280. 
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FIG. 9: Same as Fig. Elbut for the initial velocity field v = —0.1 sin(2#) 9. The shell is shown 
at times t' = 1.273, 1.280, 1.283, 1.285, 1.287, 1.2875, 1.2880 and 1.329. The earlier evolution is 
qualitatively the same as in Fig. The continuous infall of many shells with this initial velocity 



field produces the caustic shown in Fig. 11 (a 
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(b) Dark matter flows. 

FIG. 10: Cross section of the inner caustic produced by the irrotational axially symmetric velocity 
field v = -0.05 sin 29 9. 
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FIG. 12: The butterfly catastrophe has three cusps, and three points where the caustic intersects 
itself. The number of flows in each region is indicated. 
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FIG. 13: Cross sections of the inner caustics produced by the axially symmetric initial velocity field 
of Eq. {23 with gi = -0.033, and (a) c\ = 0, (b) c 2 = 0.01, (c) c 3 = 0.05, (d) c 3 = 0.1. Increasing 
the rotational component of the initial velocity field causes the tent caustic (a) to transform into 
a tricusp ring (d). 
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FIG. 14: Dark matter flows in the z = plane for the initial velocity field of Eq. (|28[) when (a) 
^ = 0,g = —0.05, and (b) £ = 0.01, g = —0.05. In (a) the particles move on radial orbits. In (b) 
the particles have angular momentum about the z axis. 
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FIG. 15: The inner caustic for the initial velocity field of Eq. (|28[) with £ = 0.005, g = —0.05: (a) 
the tent caustic in 3 dimensions; (b) a succession of constant z sections over the range < z < 0.01. 
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FIG. 16: The tent caustic of Fig. ED in (a) y = cross section and (b) x = cross section. The 
manner in which the tent pole is connected to the tent roof is described in detail in Fig. 1171 
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FIG. 17: (a)-(d) Constant z sections of the tent caustic of Fig. ^]in the region where the tent pole 
connects with the tent roof. In (a) the pole is entirely inside the roof. In (b) the cusps of the pole 
which are in the x = plane have pierced the roof. From (b) to (d) the parts of the pole and the 
roof near the y = plane traverse each other by forming hyperbolic umbilic catastrophes, one on 
the x > side and one on the x < side, (e)-(h) Constant z sections of the tent caustic near the 
hyperbolic umbilic on the x > side, for greater detail, (i) the hyperbolic umbilic at z > 0, x < 
in three dimensions. 
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FIG. 18: The first, second and third columns show respectively the z = 0, y = and x = 
cross sections of the inner caustics produced by the irrotational initial velocity field of Eq. (|28|) for 
(51,52) = (-0.033, -0.033), (-0.04,-0.02), (-0.05, 0), (-0.06, 0.02) and (- 0.067, 0.033) in five rows 
from top to bottom. During this sequence, the caustic transforms from a tent symmetric about 
the z axis (first row) to a tent symmetric about the x axis (last row). The roof of the caustic in 
the first row becomes the pole of the caustic in the last row (blue) and vice- versa (red). 
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(e) (f) 

FIG. 19: z = sections of the inner caustics produced by the initial velocity field of Eqs. (|18|) for 
ci = C2 = 0, £ = 0.01, g = —0.05, and increasing values of c 3 : (a) C3 = 0, (b) c 3 = 0.005, (c) 
c 3 = 0.01, (d) c 3 = 0.015, (e) c 3 = 0.025, (f) c 3 = 0.1. 
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FIG. 20: Inner caustics produced by the initial velocity field of Eqs. (|18|) for £ = 0.01, g = —0.05, 
ci = C2 = and increasing values of C3: (a) C3 = 0, (b) C3 = 0.015, (c) C3 = 0.03, (d) C3 = 0.12. 
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